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Analysis of the effective thermal conductivity of ceramic coatings and its relation to the 
microstructure continued. Results (obtained in Task 1) for the three-dimensional problem 
of heat conduction in a solid containing an inclusion (or, in particular, cavity - thermal 
insulator) of the ellipsoidal shape, were further advanced in the following two directions: 

- closed form expressions of H tensor have been derived for special cases of 
ellipsoidal cavity geometry: spheroid, crack-like spheroidal cavity and needle shaped 
spheroidal cavity; 

- these results for one cavity have been incorporated to contrast heat energy potential for 
a solid with many spheroidal cavities (in the approximation of non-interacting defects). 
This problem constitutes a basic building block for further analyses, since the 

ellipsoidal shape covers a variety of practically important pore geometries. 



The problem is formulated as the determination of the change in thermal conductivity 
due to inclusion. Namely: 


AQ = H G 


( 1 ) 




where AQ is the heat flux change per reference volume V , G is the far-field temperature 
gradient and second rank tensor H is a function of the inclusion shape and the inclusion 
conductivity. 

Mathematical considerations based on the analysis of the ellipsoidal shapes in the 
framework of Eshelby-type theory and on utilization of Green's function for the heat 
conduction problem in an unbounded medium, show that H has the following form: 


H = —(k* — * 0 )(A 1 // + A 2 mm + A^nn) 
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where V* = 


An 
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is the volume of the ellipsoidal inclusion with semi-axes a^, 02, <33 


with unit vectors l,m,n; k* and )cq are conductivities of the inclusion and of the matrix, 
correspondingly. Coefficients A lt A 2 , A 2 are given in terms of elliptic integrals (see Report 

for Task 1). 

In the case when the inclusion is a spheroid (aj = 02 = a). tensor H takes the form 
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or, in components: 
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where it is denoted: 

y = a]a 2 - aspect ratio of the spheroidal inclusion, 

n = n^i + n 2 e 2 + n 2 e 2 - unit vector along the axis of symmetry of spheroid, 

I = + e 2 e 2 + e 2 e 3 - unit second rank tensor, 

fo(y)= s(y) = ~lA arctan-^y 2 - 1 (for oblate shape, y > 1), 

2(1 -r) i/r-i 

2i/i-y 2 1 - V 1 - y 



(for prolate shape, y < 1). 



In the case of spheroidal cavity (insulator, k* = 0) tensor H is as follows: 
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• In the case of thin spheroidal cavity ( y » 1): 

H = -—kS -(I-nn) +— nnj 

V ^[l + 7r/(4y) v k J 

• In the limit of a circular crack : 
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In the case of needle-shaped spheroidal cavity ( y « 1): 
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In the approximation of non-interacting cavities (each cavity experiences the influence 
of the same far-field temperature gradient G unperturbed by the presence of other cavities), 
the heat energy potential AQ for a solid with many cavities is obtained as follows (in 
terms of derived tensors characterizing i - th cavity): 


AQ = -G- 
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For example, in the case of spherical cavities we have: 
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An = ~pkoG-G = -jpko(G{ + C| + C?) 
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where parameter p is the conventional porosity. 
In the case of circular cracks we have: 



( 11 ) 


AQ = -—k$GaG 


(11a) 


where a is the second rank crack density tensor (well known in problems of effective 
elastic properties of cracked media). 
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